1, 2008 5:31 WSPC/INSTRUCTION FILE charm06 



CHARMONIA FROM LATTICE QCD 

GUNNAR S. BALI* 

Institut fur Theoretische Physik, Universitdt Regensburg, Universitdtsstr. 32 
93053 Regensburg, Germany 

gunnar. balilSphysik. uni-r. de 

Recent lattice QCD results on charmonium properties are reviewed. I comment on 
molecules and hybrid states as well as on future studies of states near strong decay 
thresholds. 
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1. New States 

Recently, for the first time in over twenty years, several narrow hadronic resonances 
have been discovered. Examples are an T Z?-wave, the B Cy the r]' CJ the h c and also new 
charmed baryons. Some findings were particularly embarrassing for theorists, such 
as the D* (2317), the Z?* 1 (2460) and the AT(3872), whose masses and properties were 
not anticipated. Of course, in retrospect everyone can explain anything. The power 
of post diction is symptomatic for our lack of understanding of strong interaction 
dynamics. This is further illustrated by the now fading pentaquark-hype. What are 
the relevant degrees of freedom for the classification of hadronic excitations? Are 
angular and radial excitations of quark model states sufficient to obtain a qualitative 
picture or do we have to include hybrid- and molecules (i.e. meson-meson and/or 
tetraquark states) as well? (How) can we make quantitative predictions of masses 
and strong decay rates? 

Charmonia offer an exciting arena for improving theoretical tools and experimen- 
tal techniques. In some aspects charm quarks can be classed as "heavy" . Nonethe- 
less, with relative valence quark velocities v > 1/2 and binding energies A > m c /3, 
it is clear that corrections to an effective heavy quark field theory description can be 
large: ultimately there is no short-cut to the inclusion of non-perturbative low en- 
ergy QCD effects. What we learn from charmonia will to a large extent also apply to 
light hadron spectroscopy. However, (potential) non-relativistic QCD ((p)NRQCD) 
expectations can still serve as a starting point. Also potential models with phe- 
nomenologically fitted parameters appear to describe many features reasonably well 
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Fig. 1. The eight recently discovered charmonium states. Some can be interpreted as cc quark 
model states (standard) while others do not appear to fit into the expected pattern (?????). 



and it needs to be clarified to what extent these are connected to QCD as the fun- 
damental theory. 

Due to the spatial compactness of charmonium wavefunctions we encounter the 
luxury of finding many narrow resonances below strong decay thresholds. Also above 
such thresholds the situation is still much cleaner than it is in the light quark sector. 
While this applies even more so to bottomonia, in the T case, production in e + e~ 
machines is limited to vector states and their decay products. Not even the r)b meson 
has been detected thus far! Furthermore, the BB threshold region is elusive. Some 
information could be gained by scans above T(4S) and runs at T(10860)|T(11020) 
in a future super- £? factory. Nonetheless, the prospects in T spectroscopy are quite 
bleak due to the small hadronic production cross sections and the unavailability of 
a top-meson that could produce T states in its decay. In the charm-sector however 
many new states have been discovered recently, some in pp collisions at FNAL, 
others in -0 production and decays at CLEO and many by the B factories, mostly 
in B decays. A summary of these recent findings is presented in Fig.^ 10 resonances 
have been discovered 1974 - 1977, none in 1978 - 2001 and eight (!) since 2002. For 
recent reviews of the properties of these new states, see e.g. RefN 
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Fig. 2. The quenched charmonium spectrum (ColumbiE0, CP-PAC£0), glueballiPland spin-exotic 
cc-glue hybridJ^j overlayed with the experimental spectrum. 



2. Spectroscopy: rif — and rif > 

Most hadronic resonances undergo strong decays. This situation is hard to study on 
the lattice and requires repeating the simulation on different (and large) physical 
volumes. Fortunately, most charmonia are very narrow. The states above 3900 MeV 
however have widths ranging from O(10) MeV to O(100) MeV. A model with strong 
decays switched off can serve as a starting point for investigations of the more 
involved realistic situation. This can be achieved by switching sea quark effects off 
completely. Studies in this so-called quenched approximati on have be en performed 
by a variety of collaborations, with an anisotropy parameter _ _ _ f $ = Q« / Q * > 1 



tall 011111 -211 alf 41 



and keeping spatial and temporal lattice resolution the same 
a t . 

For charmonia, 3> M, where the mass M is calculated in units of a^ 1 

and the wavefunction extent r in units of a s . Therefore, introducing an anisotropy 
appears to allows for a smaller number of lattice points and hence a reduction in 
the computational effort. However, quantum corrections considerably reduce the 
naive gain factor. Also there is some overhead from the nonperturbative matching 
of fermionic and gluonic anisotropies, which is particularly delicate when sea q uark s 
are included. Nevertheless, a pilot study with rif — 2 and £ w 6 already exists.^^ 

We summarize the quenched situation in Fig.0 The ColumbicPand CP-PACS^ 
results are extrapolated to the continuum limit. Similar results for a smaller number 
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of states have been obtained by the Guangzhou grourJ^Hand QCD-TARO.Iil It is 
reassuring that the latter study on isotropic lattices confirms the anisotropic results. 
The scale r^ 1 w 394 MeV is set from potential models. A systematic scale error of 
about 10 % should be assumed, for stable states that do not mix! One often thinks 
of the charm quark as heavy in the sense that its mass is bigger than typical QCD 
binding energies. However, the spectrum of states entirely made out of glue covers a 
similar energy range. To illustrate this we have superimposed these glueball^ onto 
the Figure. 

Note that for the lattice results displayed, effects of diagrams with disconnected 
quark lines have been neglected. One might expect OZI violating contributions 
from these, for instance for states that lie close to glueballs with the same quantum 
numbers or close to decay thresholds. Also (tiny) effects from the axial anomaly 
or from r\ c — rj' mixing cannot be excluded. The impact of such diagrams on the 
vector and pseudoscalar channels has been studied by McNeile and Michael" 1 and 
by CP-PAC^SI and found to be insignificant within errors of about 20 MeV. 

The Figure also contains the lightest two spin-exotic cc-gluon hybridsP Again, 
the results agree with other studies j gLJSl It is also possible to employ hybrid sources 
with non-exotic quantum numbers. States prepared in this way are found t o de- 
cay rapidly in Euclidean time into the respective non-hybrid ground stateJ^^^ 
which is hardly surprising, given the fact that one expects an energy gap of about 
1.4 GeV when solving the Schrodinger equation with ground state and hybrid lat- 
tice potentials as input. However, the extent of spe ctros copic hybrid components 
in charmonium wavefunctions deserves further studyP^El Similarly, it is difficult to 
disentangle S and D waves in lattice simulations of the 1 sectorN 

When sea quarks are switched on, would-be glueballs, standard charmonia and 
DD molecules will mix (as will exotic hybrid charmonia and molecules). This has 
not yet been studied on the lattice. In addition, one would expect the fine structure 
to be particularly sensitive to the sea quark content. For instance, the S wave 
singlet-triplet splitting in NRQCD reads, AM = cp / (6™%) (ip\Vi\t(>) + •■• where 
cpVi = (327r/3)a s <5 3 (r) to leading order in perturbative QCD. Due to the short- 
distance nature of this term one would expect significant finite a effects. Moreover, 
with sea quarks included, the running of a s (q) will slow down such that at high 
momenta/short distances a s will maintain a relatively larger value. In RefM-'-'this 
effect is estimated to amount to 30-40 %. 

Indeed, setting the scale from ro, a splitting of only 77(2) (6) MeV is foundP^ 
When adjusting the lattice spacing to the IP — IS" splitting instead, the central 
value increases to A M ss 89 MeV, still short of the experimental AM « 117 MeV. 
This result is consistent with Refs^ 6,!2, while in small volume studies at one lattice 
spacin^^l and at a to o sm all charm quark masJ^ higher results have been obtained. 
The FNAL and MILO^l collaborations indeed find the J/ip — r] c splitting to increase 
to AM = 94(1), 101(2) and 107(3) MeV at lattice spacings a" 1 w 1.1,1.6 and 
2.3 GeV, respectively, in studies witn^n/ =2 + 1 sea quark flavours. Note however 
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that there are theoretical doubts about the validity of MILC's fermion approach. 
Also, disconnected quark line diagrams have not yet been included. While their 
impact is likely to be small effects of up to 20 MeV cannot be excluded!^ 



3. Radiative Transitions 

In order to match experimental resonances to theoretical predictions ultimately one 
has to go beyond the mass spectrum. For instance one would expect hybrids and 
molecules to couple differently to decay and production channels than states with 
quark model dominance. Gaining information on strong decay rates from lattice 
simulations, without additional assumptions, is complicated. However, calculating 
electromagnetic (EM) decay and transition rates is in principle straight-forward. 
For narrow states such predictions can directly be confronted with experiment. For 
broad states above threshold, where data on EM transitions is lacking, a knowl- 
edge of the wavefunction will reveal qualitative information on possible production 
mechanisms. 

QCD-TARCP] calculated electric quarkonium wavefunctions. In the non- 
relativistic approximation, the EM decay constant is to leading order propor- 
tional to the wavefunction at the origin. Dudek and Edwards directly calculated 
two-photon decay rateJ^ as well as the decay constants^ /jm, f^>>, /jj c and f v > r . 
These studies have been performed at one lattice spacing on small volumes in 
the quenched approximation and with a charm quark mass about 5 % smaller 
than the physical one. In particular the calculations of this pilot study of EM 
transition rates are extremely promising. For instance the Ml transition rate 
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T(J/ip — > 7y c 7) cx q 3 /(™j/psi + m Vc ) 2 a {s \V(0)\ 2 contains the term V(Q 2 ) at Q = 0. 
This decomposes into a kinematical factor, multiplied by the three point function, 
(Vc(p') |j M (0) I J/^ip)) 7 that can be calculated on the lattice. In Fig.|3]we display the 
lattice prediction-^ of V, obtained for different momentum transfers Q 2 = (p' — p) 2 , 
together with the Q = value derived from the experimentally measured transi- 
tion rate. This expectation was phenomenologically corrected for the mis-adjusted 
lattice charm quark mass (lat. mass) to facilitate the comparison: the lattice Q 2 - 
extrapolation appears very reasonable. 

The transition Xd J I^Pl was studied too as were Xco ~ ► J 1^1 an d r(/i c — > 
j] c -f) w 600 keV. The latter rate is a prediction since CLEO only detected the 
cascade tp' — > itoh c , h c —* r/cj where both individual rates are unknown. 



4. Mixing and Strong Decays 

It is evident that states close to decay thresholds like the X(3872) or the F(3940) 
(see Fig. ^ will couple to molecular DD or D S D S components, respectively. At 
present there exist no meaningful direct lattice studies of such cc <-> cqqc mixing 
phenomena. Ideally one would create and destroy states using an extended basis of 
operators, including molecules and hybrids, to study their mixing. 




In the absence of such an investigation, one might benefit from the non- 
relativistic approximation. At least qualitatively charmonia can be described by 
potential models. The QQ sector can be extended to include W potentials. Q de- 
notes a static (he avy) quark and T> — Qq. This situation has recently been studied 
on the lattice*^"" albeit only at one lattice spacing a w 0.083 fm and with n/ = 2 
mass-degenerate sea quark flavours, slightly lighter than the strange quark. We 
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denote the physical eigenstates by \n). These can be decomposed, 

\l) = cos 6 \QQ) + sin d\VV) + --- , (1) 
|2) = -sm9\QQ) +cos6\VV) + • •■ . (2) 

The resulting energy levels and mixing angle 9 are displayed in Fig. 01 as a function 
of the QQ-separation r. 

At r > r c w 1.25 fm the QQ string breaks and the ground state is dominated 
by the T>V configuration. Interestingly, also at r <C r c there is a significant four- 
quark admixture in the ground state. This 4 % probability is consistent with the 
large four-quark components found in phenomenological models of charmonia and 
bottomonia!^] In Refill the action and energy density distributions around the 
static sources have been investigated: string breaking appears to be an instantaneous 
process with no spatial localization of the light qq pair that is created. This suggests 
an 2 x 2 Hamiltonian with a transition term (E2 — Ex) sin 9 cos 9 that only depends 
on the distance r. _____ 

One can include T>T> effects as quantum mechanical perturbations ,____- starting 
from an unperturbed QQ wave function ipo. By summing over all possible interme- 
diate states i where ipi can be in both sectors, T>T> and QQ, one ends up with 
complex Hamiltonian for each state. The solution of the corresponding Schrodinger 
problem yields both: the respective energy levels and decay widths. 

A study of the sea quark mass dependence of the static energy levels is as 
yet lacking. The bands in the left plot of Fig. 0] represent an n/ = 2 + 1 spec- 
ulation. In this case string breaking will occur at a somewhat smaller distance, 
7'c = 1.13(10) (10) fm, the gap E2(r c ) — Ei(r c ) will be larger and moreover there will 
be a second string breaking threshold QQ <-> T> S T) S . It would be interesting to see, 
once the sea quark mass dependence has been established, if properties of X(3872) 
and Y(3940) can be reproduced or if the present overpopulation of the 1" sector 
can be explained by coupled channel models, based on such lattice potentials. 

5. Summary 

We have arrived at a good understanding of the spectrum of states in the mass 
region, 2.5 - 5 GeV, in the quenched approximation to QCD, including glueballs 
and cc-glue hybrid states, besides standard charmonia. However, information on 
four-quark states is still lacking and there is some uncertainty regarding the size 
of contributions due to diagrams with disconnected quark lines. Electromagnetic 
decay and transition rates have been calculated successfully in a pilot study. 

Studies of n/ = QCD set the stage for more realistic nf > simulations. 

1 

Promising rif = 2 + 1 results have been obtained, indicating a convergence of low 
lying states towards experiment, once sea quarks are included. Over the next few 
years such simulations will be extended to include more realistic sea quarks and a 
broader basis of operators to create states, including four-quark and hybrid opera- 
tors. In the meantime one has to resort to coupled channel models to understand 
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threshold effects. These can be rehned by using input from lattice simulations. Ul- 
timately, one would wish to arrive at QCD predictions of strong decay rates. This 
is very hard to achieve by a direct lattice computation but again, the predictive 
power of models can benefit from dedicated lattice tests and comparisons. Last but 
not least I hope that these exciting times will continue with further experimental 
discoveries. 
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